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SIAMHY 


V region thickness , c a 


Recent use of alloying techniques for rear contact 
fonatlon has yielded a nev kind of silicon solar cell, 
the back Surface field (BSP) cell, with abnozmlly high 
open circuit voltage (V^) and isproved radiation re- 
sistance. Several analytical models for based on 
the reverse sat uration current are formulated to ex- 
plain these observations. The zero SHY case of the 
conventional cell model, the drift field model, and the 
l 0 *-hdg£ Junction (LRT) model can predict the experi- 
mental trends, the IHJ model applies the theory of 
the lov^high Junction and is considered to reflect a 
more realistic riev of cell fabrication. This model 
can predict the experimental trends observed for BSP 
cells. Detailed descriptions and derivations for the 
models are included. The correspondences between them 
are discussed. This modeling suggests that the tinning 
of minority carrier diffusion length measured in BSP 
cells be reexamined. 

TABUS OP SHOOTS 

minority carrier diffusion coefficient, 
cs?/aec 

electric field strength, volt/cm 

minority carrier generation rate, sec*’ 1 ' 

short circuit current density, amp/ca? 

diode saturation current density, amp/an 2 

minority carrier current density, amp/ as? 

dsnege constant for 1 MeV electrons, 
particle" 1 

Boltzman constant, 1.38X10* 2 3 Joules/K 
minority carrier diffusion length, cm 
minority carrier concentration, an“3 
excess electron or hole concentration, oa*3 
intrinsic carrier concentration, cn*^ 
donor concentration in n-type material, cm* 3 
N A+ acceptor concentration in p, p+ material, cm* 3 



x arbitrary distance into the cell, cm 

Xj depth of metallurgical Junction, aa 

X subregion boundary plane location, cm 

u minority carrier mobility, ca 2 /volt sec 

t minority carrier lifetime , sec 

* flue ace of 1 MeV electrons , cm* 2 

♦ equilibria barrier potential of the low-high 
junction, volts 

gqpacripta 

p associated jlth p-type region or with 

minority carrier electrons 

p+ associated with p-*-type region 

a associated with n-type region or with 

minority carrier holes 

n+ associated with n+-type region 

IHTHODUCTIOM 

Recent efforts to improve the performance of sili- 
con solar cells have resulted in an increased open 
circuit voltage (Y^) for 10 oim-cm silicon material. 
For years the V oc of the 10 oha-cm silicon cell was 
fixed at about 0.55 volts. However the use of an 
alloying fabrication technique (1) for rear contact 
fbmtlon has resulted in the bach- surface field (SSf ) 
cell having voltages as high as 0.58 volts. In addi- 
tion, imp ro v ements in short circuit current, efficiency 
and fill factor have also been observed. Unlike the 
conventional 1C oho- cm cell, the of the BSF cell 

appears to be independent of cell thickness at least 
over the thickness range of 100 to 400 micrometers. 
These improvements In cell performance are directly 
attributed to the presence of a heavily doped p-type 
layer (referred to as the p+ layer) formed during rear 
contact fabrication. Measurements of the effective 
diffusion length of these cells using various techni- 
ques have produced uncertain Information. In some cases 
the effective diffusion lengths are in excess of values 
measured for conventional cells. 


q electronic charge, 1.6*10* ^ coulombs 

R minority carrier recombination rate, sec* 1 

b surface recombination velocity, cm/sec 

S normalized surface recombination velocity for 

low-high junction model, dimensionless 

t total solar cell thickness, cm 

T temperature , kelvin 

v open circuit voltage, volts 

oc 


Analytical modeling for the BSF cell has not re- 
ceived as much attention as the experimental work. 

Thus the purpose of this paper is to present calcula- 
tions of v oc based on conventional, drift field, and 
low-high junction one dimensional solar cell models and 
to test these calculated results against the experi- 
mental data. The model for the conventional cell as- 
sumes an arbitrary front and back surface recombina- 
tion velocity (SRV). This model has four limiting 
cases: 1) the infinite 8HV case, 2) the zero SRV case, 

3) thr SRV equals diffusion velocity case, and 4) the 
infinitely thick cell. The drift field model for V oc 
is based on the same assumptions and boundary condi- 
tions as Rblf’s drift field model used for the calcula- 
tion of collection efficiency (2). The new model to be 
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proposed here is the low-high Junction (LHJ) model* 

The requirement for a IH J rear contact structure (p*p+) 
In the BSF cell results In boundary conditions that 
reflect the minority carrier blocking proper. .. of the 
LOT first reported by J. B. Gunn (3)* The UU m lei 
presented here differs from that derived by button and 
Whittier (U) by neglecting space charge recombination 
effects. However this work extends the work of ref. h 
by including the p+ region geometry fhctor. 

Although LHJ effects are well-known and have been 
applied to the modeling and analysis of semiconductor 
rectifiers (5), this paper represents a first attmqvt 
at applying this model to photovoltaic energy conver- 
ters. 


The approach taken in the modeling is to relate 
the V to the diode saturation current (Iq) and 
short circuit current (l gc ). Expressions for Iq are 
derived assmlng boundary conditions for low Injection 
minority carrier transport with boundary conditions 
appropriate to each model. The performance of all the 
models and cases Is cospsred with the experimental 
Yoc versus thickness data. Finally, the radiation 
degradation behavior of Voc i» calculated using the 
IBJ model. 


If it is asswed that the current transport in a 
solar cell is dominated by min o rity carrier diffusion 
processes, then the open circuit voltage (Vqq) can be 
related to the short circuit current (Ig C ) and the 
diode saturation current (Iq) by the following simple 
expression 



Hence an analysis of the Voc can be made far the con- 
ventional, drift field, and low-high Junction models 
of the solar cell using expressions for Ip which re- 
flect the cell structure and properties unique to each 
model. These expressions are given below. Calcula- 
tions were made using a high speed digital counter 
and a programmable desk top calculator . All calcula- 
tions were made for short circuit density of about 
4o ma/<=? and a diffused region saturation current 
(Ioq) equal to 1.2XL0" 11 * amp/ on 2 . 

Conventional I Q Expressions 


A caamon view of the conventional n on p solar 
cell structure U one in which the p-baae dopant and 
properties ere uniform, the device dimensions are fi- 
nite , and contact effects are described by the surface 
recombination velocity (SRV). The Iq for this gen- 
eral case has been derived by Me Ke Ivey (6). This ex- 
pression is given by equat^n (2a) for the n and p 
regions. 
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The SRV is represented by s v and L is the bulk dif- 
fusion length. The other parameters have their usual 
meaning. Bote that the product of the diffusion length 
(la or I*) with the reciprocal of the bracketed terns 
defines an "effective' 4 device diffusion length. This 
equation however reduces to the following approximate 
forms. 


l) Infinite SRV Case . When s r and 6f are 
very large as would be the case for ohmic contacts, 
equation (2a) reduces to tte familiar hyperbolic cotan- 
gent form. 




qn 2 D 


(2b) 


This case gives a reasonable description of the ~ase 
region of the conventional cell structure . 

2) Zero SRV Cane . If the SRV’s are assured to 
be zero, (2a) reduces to a second form now containing 
the hyperbolic tangent geometry factor. 


(2c) 


This case lilies oon-otmic contacts and does out re- 
flect the practical solar cell. 

3) Diffusion Velocity Case . If the SRV for the 
front and rear contacts are equal to the minority 
carrier diffusion velocities (D/L) in each base region, 
the bracketed terms of equation (2a) reduce to unity. 
This approximate Iq is identical to the infinite base 
case given below. 

4) Infinite Base Case , When the base region 
thicknesses are such greater than the bulk diffusion 
lengths, the bracketed terms in equation (2a) again 
reduce to unity. The resultant Iq given by equation 
(2d) is the well-known expression derived by Shockley 
(7). 


Vb 


(2d) 


This case and the preceding case predict the I Q to be 
lndepeixlent of base region thickness. 


Drift Field Model Ip 

Practical solar cells are fabricated by diffusion 
and alloying processes wnich usually result in regions 
having significant nen-uniform Impurity distributions. 
It Is easily shown that such distributions give rise to 
an electric field which can significantly influence 
minority carrier transport in those regions. Kleinraan 
(8) and Wolf (2) have calculated the cell I sc using 
an expression for the collection efficiency which con- 
tains drift field effects. The cell was mathematically 
divided into subregions which can be assigned separate 
material properties and dimensions. For example Wolf 
used a four layer model in which exponential Impurity 
distributions could be assigned to each regic.1. The 
Important features of this model are that the minority 
carrier density and the current density are continuous 
at the subregion boundary. The I 0 for this structure 
has been derived (9) and includes not only finite geo- 
metry and contact effects , but also the effect of sub- 
region drift fields. The p-base component of the n on 
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p solar cell containing a drift field only In a thin 
region In front of the rear contact Is given by equa- 
tion (3). The p-type base subregion nearest the Junc- 
tion is uniformly doped and therefore field free. 
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The parameter S is defined as a normalized surface 
recombination velocity. It contains HIT Carrier 
factor and the nobility and i if fusion 1-r.gt.*. 

ratios for tfte p end p+ region. S else contains a 
geometry factor for the p* region, the SRV for the 
aetal-p+ contact (s r ), and the p+ region xhickness 
(Wp+). The total Tq must include the n-base compo- 
nent I~. In the simplest case Io- i3 described by 
the drift field model given by equation (2). riowever 
for the range of variables considered here, I r Is 
always significantly greater than for lir ohm-ca 

material. Therefore Iq„ is neglected. However the 
importance of Iq and Its effect on the performance 
of low resistivity and deep diffused cells is recog- 
nized. 


A similar expression Is written far the n-type diffu- 
sed region with the appropriate change in subscripts 
and by assuming an erfc donor distribution. The drift 
field region Is nearest the junction for this case. 

It can he seen that equation (3) will reduce to the 
general SRV case given by equation (2a) if a single 
layer, field free model is assumed. 

Low- High Junction (Iflj) Model Iq 

For the UU model, the subregion boundary, unlike 
the drift field model, represents an abrupt transition 
between the thick, lightly-doped p- region (i.e. , sub- 
region nearest the Junction) t»nd the thin heavily- 
doped p+ region. Both regions are assumed uniform in 
dopant and materials properties. It has been shown by 
Gunn that this Junction has the unique property of 
preventing excess minority carriers in the p-base from 
entering the p* region. Majority carriers pass freely 
across this junction. The minority carrier concentra- 
tions at the space charge region edges depend on the 
barrier potential (v) associated with the low-high 
junction as shown In equation (Ua ) 

", - V eXP S (4a) 


In addition, the diffusion currents at the p+ ani p 
region space charge edges are proportional to the gra- 
dients of the minority carrier densities there, and 
given by equation (4b) 


qD 7 n 
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(4b) 


Equation ( Ua ) represents a subregion boundary condi- 
tion which differs significantly from that assumed in 
the drift field model. Pecombinat ior effects and 
thickness of the L&J space charge region arc assessed 
negligible. The I 0 using this model is derived in 
detail in the Appendix. The result for the p-side la 
given by equation (5) where 1 is the cell total 
thickness 


REVIEW OF EXPERIMENTAL DATA 

Published measurements on the high voltage (BSF) 

10 ohm-cm solar cell are assembled ..ere to provide a 
comparison with the theoretical predictions. The data, 
shown in figure 1, are grouped into t-*o categories as 
follows: dependence of V QC on cell thickness (1) 

and the response of v oc to 1 MeV electron exposure 
(1, 10). The data of ref. 1 for tb e high voltage cell 
is shown in fig. la. Voltages as high as 0.59 v are 
measured and are approximately independent of cell 
thickness in the 100 to 300 micron range. Vqc varia- 
tions of ± 5 mv at each thickness are not unreasonable - 
More recent unpublished data suggest tnat constant vol- 
tages as high as 0.60 volts are achievable. These data 
are in striking contrast to these shown for the conven- 
tional cell, where thinner cells have lower voltages. 
Apparently the only difference in the fabrication be- 
tween the high voltage and conventional 10 ohm- cm cells 
is the additional alloying and/or diffusion and sinter- 
ing associated with the rear contact of the high vol- 
tage cells. The radiation behavior (l) of various 
thicknesses of high voltage cell is shown in fig. lb. 
Also shown are data for conventional cells of ISO and 
300 micron thickness. The initial voltages are higher 
for the high volta-e cell compared to 4 he conventional 
cell of comparable thickness. In addition the high 
voltage cells maintain a higher voltage for larger 
flue nee 8 than do the conventional cells. Furthermore 
the thin high voltage cell maintains a superiority over 
the thick high voltage cell. At high fluences the mag- 
nitude of the voltage and the radiation degradation 
rate of all high voltage cells are the sajr.e as fir t r.e 
conventional cell. Also shown in the figure are un- 
published data (10) for a 100 um high voltage cell 
measured by JPL. 


R££ JLTS 


A st^ly of the open circuit voltage d*pend»r.^e 

on cell structure has been mede using equation (1) and 
equations (2), (3), and (5). 
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Conventional Model 

Figure 2 shows the dependence of V 0C on p-base 
diffusion length for the infinite base case where +he 
requirement is that Wp »L^. As can be seen voltages 
in excess or 0.59° volts are predicted providing the 
p-base diffusion length is greater than 1000 micro- 
meters. This cell would be at least 0.3 cm thick in 
or er to obey the cell thickness/ infusion length re- 
quirement. Because the model assumes an infinite 
(> 0*3 cm) thickness at all diffusion lengths, there is 
of course no v oc dependence on thickness. 

Voltages calculated for a range of cell thick- 
nesses and p-base diffusion lengths using the finite 
width model are shown In fig. 3. The lotted, solid, 
and dashed lines represent calculat ions for tte infi- 
nite SKV, the zero SRV, and the special ca*e where the 
SHV ir equal to the diffusion velocity (C/L), respec- 
tively. Voltages that are independent of cell thick- 
ness are predicted only for the latter case. For very 
thick cells and a specified diffusion length, all case? 
predict voltages that ore approaching the limiting 
values of the Infinite width case. This is not unex- 
pected since the further away the rear contact is from 
the junction the less it will lnfl^nce the minority 
carrier distribution in the p-base. As tne cell be- 
comes thinner , the predicted voltage trends become 
very different. 11 la seen that the zero SKV case 
predicts voltages for thin cells that are well in ex- 
cess of those predicted by the ether two cases. In 
addition only the zero SRV case will predict voltages 
above 0.50 volts provided the p-base diffusion length 
is substantially greater than 200 pm. 

Drif t Field Model 

The results of the voltage calculations for the 
drift field model are shown in figs, 4 and 5. A 0.5 
pm thick p+ region in which an exponential acceptor 
concentration profile exists was assumed. For the re- 
sults shown in fig. 4, the p+ region diffusion length 
was assumed to be 1.0 um. It can be seen that the 
magnitude of the voltage depends on cell thickness and 
p-base diffusion length, deleft* 0.54 volt (i.e. , 
approximately 200 micrometer diffusion length) the 
Voc i» nearly independent of cell thickness. Above 
0.54 volts the voltage Increases with increasing cell 
thickness similar to the infinite SRV case. 

Figure 5 shows the variations of the voltage for 
values of p+ region diffusion length between 1 and 
about 220 pm. Calculations were made for p region 
diffusion lengths of 50, 200, and 600 pm. The p+ re- 
gion width is again 0.5 pm. As shown, the diffusion 
lengths associated with both regions control the mag- 
nitude of the voltage and its dependence on thickness, 
large values of p+ region diffusion lengths are neces- 
sary to produce high voltage 3 in thin cells. The vol- 
tages predicted for cells where the p-base diffusion 
length is much smaller than the base thickness are 
Identical to those shown in fig. 2, the infinite base 
case. 

Low- High Junction Model 

The results of voltage calculations using the low- 
high junction model are shown in fig. 6 for a range of 
p-base diffusion lengths. In practice, there is a 20 
mv spread in voltage reflecting the actual 7 to 14 
olm-<2& resistivity variation for nominal 10 oils- cm 
material. Silicon in this resistivity rang'* is called 
n ominal lo ohm- cm. A uniform p+ region acceptor con- 
centration of 1.3Xip*9/em3 was assumed yielding a 
W ratio of 10“**, The p+ region thickness and 


diffusion length are 0.5 -m and 1.0 ^ respectively. 
Contrary to the drift field model, increases in tne 
p+ region diffusion length up to 200 jr had no in- 
fluence on the V oc , The voltages for thick cells as 
asymptotically approach those cf the infinite base case 
since W/L»l. For the conditions assumed, t:.e snaps 
and magnitude of the V oc /thickness curves are identi- 
cal to those calculated for the zero SRV case. I. us, 
tne zero 5K7 case, can be used, as an approximat icr. cc 
the LHJ model, to calculate cell v 0c with p-region 
diffusion length as the only variable . As in t r e pre- 
vious models, the magnitude of the voltage is deter- 
mined primarily by the p-base diffusion length v-ich in 
itself yields limited information. Ir the LKf r.oiel 
the parameter S, which incorporates properties of the 
LHJ and the p regions, assumes major importance in 
determining cell performance. 


Figure 7 shows a plot of the voltage against tht 
magnitude of the normalized surface recombinat icr. velo- 
city S defired by equation (5b). The calculations 
are made for several values of cell thickness for a 
p-base diffusion length of 574 urn. This figure shows 
clearly the influence of £ on the magnitude of the 
voltage. For reference, the curve of figure 6 corres- 
ponding to the 600 um diffusion length has an S equal 
to about 10'^. For small values of S (less than 
unity) , the voltage asymptotically approaches the val- 
ues predicted L* the zero Spv case shown in fig. 3 • 

For large values of S, the v oc approaches that pre- 
dicted bv the Infinite SKV case of fig. 2. Thus fig. 7 
shows graphically the manner in which the LHJ model re- 
duces to both the zero SRV and infinite SRV models. 

For intermediate values of S in the range cf 0.1 lo 
10, the voltage varies greatly between the two limiting 
cases. For S equal to unity, the voltage is inde- 
pendent of cell thickness and equal in magnitude to 
that predicted by the infinite width case . 


Figure 8 shows the calculated voltage response to 
1 MeV electron exposure . Response is calculated using 
the LHJ model for cells of various thickness. Calcula- 
tions for the conventional cell of 100 and 4 00 ua 
thickness using the infinite SRV case is also shown. 

An initial p-base diffusion length of 400 um is assumed 
for all cells. The voltage degradation is simulated by 
assessing that radiation damage alters the diffusion 
length through lifetime degradation. The diffusion 
length in the Iq equations is replaced by the follow- 
ing expression containing the radiation damage coeffi- 
cient (K) and the initial diffusion length (Lp(o)): 


1 

i*(*) 


- 5 *— ♦ K 0 

L p«» 


( 6 ) 


$ is the electron flnence arid K is equal to 9 ^10“^ 
(electrons" 1 ) (11) for 10 ohm- an material. At zero 
flue nee the thin high voltage cells have higher voltages 
compared to the thick cells. This superiority is main- 
tained to about 1X10 1 5 e/enr. Beyond this all cells 
degrade at the same rate as the conventional cells. 


DISCUSSION 


Based on the above data, the infinite SRV case 
can explain the experimental trends of the conventional 
cell. As will be discussed below, the trends of the 
BSF cell are predicted using either the zero SRV case, 
the drift field model, or the IiU model, iiowever in 
all casts the "goodness of fit" is critically dependent 
on the nagnit’’ 3e of the bulk diffusion length use 1 '* 
the calculations. These diffusion lengths are not in 
agreement with those determined by experiment. Some 
reasons for this will be discussed. 
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Caspar laon of fig. la and fig. 5 shows that tte 
infinite SHV case predicts the correct V oc thickness 
trend for conventional cells. However the magnitude 
of the voltages for cells of different thicknesses can- 
not be predicted using a single value for the p-base 
diffusion length. The voltages for thin ceMs are pre- 
dicted only for a diffusion length that is lower than 
for thick cells. For example , from fig. la, conven- 
tional cell voltages of 0.55 and 0.5J V are obtained 
for cell thicknesses of 300 and 100 ja respectively. 
Calculations using the infinite SHV model require bulk 
diffusion lengths of 300 and <50 am respectively for 
these two cases. Cell measurements also Indicate a 
reduction in diffusion length. The mechanism of this 
apparent reduction in diffusion length is not under- 
stood. It is possible that reductions In diffusion 
length are effected through impurity introduction dur- 
ing processing or by cell geometry changes. 

Of more present interest are the voltages observed 
In the BSF cells. It is not possible to predict either 
the magnitude of the voltage or its trend with thick- 
ness for BSF cells using the Infinite SHV case. The 
diffusion velocity case yields voltages above 0.58 V 
only if the p-base diffusion length is greater than 
600 m- In this case the voltages are independent of 
cell thickness. Also, the zero SHV case yields vol- 
tages above 0.58 V for diffusion lengths above 300 \m 
depending on cell thickness. The thicker the cell, the 
greater the necessary diffusion length. However these 
cases provide li-tle information which is traceable to 
fabrication ~zz3 zes tne influence of cell 

structure properties on cell performance. Because of 
this, these cases will not be considered further. 

The drift field model and the LHJ model also pre- 
dict high voltages provided the V/L ratio of the p-base 
Is less than unity. It should be noted that the drift 
field model approaches the voltage/thickness trend of 
the zero SHV case of the LHJ case only for high p + 
region diffusion lengths. In addition, these models 
both contain parameters relating to controllable cell 
fabrication features, such as p+ region width and do- 
pant concentrations, as well as carrier mobility and 
lifetime. As shown in fig*. 5 and 7, these features 
can significantly influence calculated cell performance. 
Cases of observed independence of cell voltage on 
thickness are predicted by both these models. For the 
drift field model, the voltages will be nearly indepen- 
dent of cell thickness only for a particular value of 
p + -region L which is about 2.5 urn for the O.58 V 
case. The IHJ model can predict the experimental 
V 00 /thickne88 treads by applying one of the following 
approximations to equation (5): 

1 . *p»h> 

2 . 8*1 

3. Diffusion length dependent on cell thickness. 

The first approximation is trivial and yields the in- 
finite base case. The BSF cells do nut physically con- 
form to this case. The second approximation, however, 
imposes specific relationships between certain cell 
properties. Equation (5b) for S can be reduced to 
equation (6) by assuming an ohmic contact to the p+ 
region 


( 6 ) 


For S b l, an appropriate relationship must be estab- 
lished between the p + region geometry, the diffusion 
velocities (D/L), and the IHJ barrier height (Na/H a +). 


It le easily seen that there is no unique set of para- 
meters which satisfies the condition. Therefore the 
S * 1 approximation would seem to necessitate a dif- 
ficult interaction and/or counterbalancing of multiple 
cell parameters. 

Examination of fig. 6 suggests that if a line 
corresponding to C .58 volts, for example, were ir avr. 
parallel to the thickness axis, a ur.icje value :: :::• 
fusion length would correspond to each cell thickness. 
Using this method, a 150 jd thick cell would rea.iire a 
300 urn diffusion length while a 500 -m thick cell would 
require a 5 00 um diffusion length. 

This trend of decreasing diffusion length with 
decreasing thickness is also observed in conventional 
cells as shown previously. 1 if fusion length mess -re - 
ments of lies (12) on wafers ranging from Loo r 
several centimeters thick also suggest that diffusion 
length does indeed decrease monotonicaliy as sample 
thickness decreases. The measured L varied from *30 
um for the LOO um thick wafers to over 1000 um for the 
very thick sanples. All samples had oimic contacts. 

Ho data presently exists for samples having "zero SRV”- 
like contacts as might be expected from a high voltage 
(BSF) cell rear contact. Therefore, from the trends 
seen in cells, a decrease uf diffusion length vior. de- 
creasing cell thickness appears to be reasonable, al- 
though no mechanism has been advanced as yet to ration- 
alize the effect. 

The requirement that high voltage cells have dif- 
fusion lengths that are substantially larger than the 
200 um measured for conventional cells also is of con- 
cern. However, in addition to the experimental evi- 
dence of lies (12) that long diffusion lengths exist, 
other theoretical and experimental bases can be found 
in the literature. The published work by Zinmeraan (13) 
suggests that the IHJ causes an excess minority carrier 
distribution in the p-base that differs significantly 
from the simple exponential profile that occurs in a 
conventional cell geometry. Therefore, it is suggested 
that methods of measuring diffusion length that rely on 
semi- infinite cell dimensions and/or o!roic contacts, 
although used on the LHJ structure, cannot be in the 
usual straightforward manner to yield the bulk diffu- 
sion length when a IAJ is present. In addition there 
is soots experimental evidence from the work on epi- 
taxial structures (1L) that the presence of a LHJ 
structure results in a higher measured device diffusion 
length compared to similar devices without the n+-n or 
p+-p contact. This work suggests again that the usual 
methods of determining device diffusion length must be 
critically reviewed in light of the carrier accumula- 
tion (or blocking) property of the LHJ. Furthermore, 
the effect of finite sample dimensions, particularly as 
regards the influence of surface effects on the bulk 
diffusion length must also be considered when inter- 
preting a diffusion length measurement. 

Let us now turn our attention to other features 
of the LHJ model which are useful for interpreting the 
experimental results. The LHJ model appears to have 
vide applicability and permits a degree of insight into 
the influence of fabrication methods on cell perfor- 
mance. For example, the LHJ model has been applied to 
the pre- and post -irradiation performance of the epi- 
taxial solar cell (15) with considerable success. Ir. 
addition, the predicted radiation performance of tnick 
and thin PSF cells, shown in fig. 8, are also in rea- 
sonable agreement with measured trends. 

It is known that not all cells fabricated using 
the BSF process exhibit the high performance. The LHJ 
model via the parameter S (equation (5b)) suggests 
that in some cases the low voltages may be traced to 
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an imperfect Ln barrier (NaZ^A*) and unoesiraole p + 
region properties (width and diffusion length) both 
of which can result in large S values. The parame- 
ter S can be interpreted as a dimensionless quantity 
related to leakage of minority carriers in the p-:e- 
gion across the XH barrier into the p* side. The 
greater this leakage is , the larger the magnitude of 
S. This leakage is influenced by at least the follow- 
ing three effects: 1) a poor quality UU due tc a 

diffusion anomaly, 2) an undesirable diffusion velo- 
city (D/L) ratio, and 3) a small p + region W/L ratio. 
The first effect may arise from crystallite imperfec- 
tions which enhance metallic precipitation and impur- 
ity spiking. In addition the alloying or diffusion 
process may be critically dependent on crystalline 
perfection and the presence of imparities which in 
turn would influence the p* region thickness and the 
carrier lifetime. Small values of the p* region W/L 
ratio can cause substantial increases in S even 
though the Na/Na+ ratio may be desirably small and 
the p- region diffusion length relatively large. The 
consequence of these effects is that the cell would 
behave as if it had an ohmic rear contact, and hence 
yield voltages approaching those of the conventional 
cell. 

In addition to the quality, the uniformity of the 
loir-high 1 action rust also Influence cell performance. 
The effect of localised regions of the IiU having 
widely different H^/N a + ratios, for example , is not 
known. This is true not only for the influence on 
v oc but also on the measured cell diffusion length. 

SUMMARY AND CONCLUSIONS 

The open circuit voltage as a function of cell 
thickness has been calculated using the diode satura- 
tion current expressions appropriate to several solar 
cell structure models. Expressions were derived for 
the diode saturation currents appropriate to the drift 
field model and a new model proposed here , the low- 
high junction model. Also investigated were the infi- 
nite and finite base width cases of conventional cell 
model. In the latter, the zero and Infinite surface 
recombination velocity approximations were considered. 
Comparisons were made with experimental data for back 
surface field (BSF) and conventional cells. The low- 
high Junction model was compared to tke experimental 
data for an un irradiated cell and cells irradiated 
with 1 MeV electrons. 

Based on these analyses the following conclusions 
are made: 


the infinite SftV case applied to the conv .tional 
solar cell. The requirement of a hign diffusic- lergth 
poses an unresolved dilemma which suggest a further 
exploration of the factors that influence diffusion 
length. 

5. The LHJ has been applied to the epitaxial 
solar cell and successfully predicts the raeasur - i 
trends. 

APPENDIX 

DERIVATION OF I 0 FOR THE LCK-HIGH JUHCTJC:: 

MODEL OF THE BSF SOLAR FELL 

The cross section of a solar cell with a 1-.- .igh 
junction as the rear contact is shown in fig. 9- The 
derivation of the saturation currei. for this Czli 
structure begins with the following equations basic 
to any semiconductor device analysis: 


Continuity Equations 


dn 

dt “ 

G -R + i v * Jn 
n n q 

electrons 

(1A) 

dt 

G -R - i 7 • Jp 
P P <1 

holes 

(2A) 

Current Transport Equations 


J n 

*= qp Q nE + qD Q vk 

electrons 

(3A) 

J 

P 

= q^pE • 

holes 

(UA) 


Bo is son Equation 

V • E - (p-n+Hp-up (5A) 

where 

B e - v ♦ (6A) 

Equations (1A) to (4A) may be combined following 
van Roosbroeck (16) to give the ambipolar continuity 
equation assuming space charge neutrality ( i.e . , 

* £p) 

dt 1 ■ G n‘ R n + ** Ev + D* ^ (tf») (7A) 

where 


1. The zero 8KV case of the conventional cell 
model, the drift field model, and the IiU model can 
predict the high voltages measured in BSF cells. High 
voltages are predicted only when the W/L ratio of the 
p-base is lees than unity However only the drift 
field and LHJ models describe cell performance in terns 
of parameters that are traceable to the fabrication 
process. 

2. Based on the results of calculations, and the 
present lack of experimental data , no clear choice can 
be made as to which model is more appropriate for the 
high voltage BSF cell. 

3. Calculated voltage degradation due to 1 MeV 
electron bombardment was compared with a limited amount 
of experimental data. The experimented trends of vol- 
tage degradation with cell thickness was predicted. 

h. The p-base diffusion lengths required in the 
drift field and IHJ models of the high voltage cell 
are about 2 to 3 times higher than those required in 


ambipolar mobility 


P/^n + 


ambipolar diffusivity 


P/M n ♦ n/u p 


For extrinsic material p £ p R and D 3 D n and 
equation (7A) reduces to the more familiar expression 
describing continuity of minority carrier transport 

£ • V R n * »„ E <•»> * ”, ^ <*> 

A more general statement of (iOA) can be written as- 
suming that not only the carrier density but alsu the 
mobility end the impurity concentrations are spatially 
dependent. Using equation (1A) to (5A) again, this 
equation has the form, for no external carrier genera- 
tion and a steady state mode of operation 


■$7a 
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It should be recognized that In t hie case an effective 
electric field 5 la defined by 


flpfx) ■ n 


£ 


kT 




■ E 


and an effective diffusion length L by 


. X-x - . . X-x 

cosh -j— + S 8 lnh 


cosh + S sinh — 
? P 


(iu; 


where S la defined aa a normalized recaab 4 nation 
velocity and is given by the following expression 


qB^ 
“ •*> 


•f V 


as.i 

kT .2 
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In the field free case assuming a ncmogeneous material 
equation (HA) reduces to the simple form for low in- 
jection 

J n . t*- » 0 (12A) 

n 

This equation has the general solution 

n(x) - A ocsh x + 0 sinh x (13A) 

In applying this solution to the p-baae solar cell 
containing a high low j* met ion the following boundary 
conditions are employed 

A. x » Xj, n » n^fexp 0 V - 1) — junction edge 

B. x * X, ruCX) * ii + (X)ejqp 0 ^ — subregion 

p v boundary 

flPp - <vv (x) 

C. x . W, qD + (W) - -q s n + (W) — meial - p+ 

p p p interface 

Boundary condition B is of prime importance since It 
contains the potential barrier associated with 

the high-low junction and expresses the minority car- 
rier concent rat ions at the space charge edges. The 
condition for current transport across the x * X boun- 
dary assumes that space charge recombination is negli- 
gible and that current is primarily by diffusion pro- 
cesses. This has been implied by equation (12A). Any 
excess minority carriers appearing at the p+ region 
space charge edge as a result of the biasing of the 
n + p junction are assumed to diffuse to the metal rear 
contact. The recombination of tbs excess carriers 
there is controlled by the surface recombination velo- 
city f s. 

The solution for the excess minority carrier con- 
centration for the one dimensional case is derived by 
dividing the p-region Into two subregions the p-base 
and the p 4 base. The subregion boundary occurs at 
x * X. Both regions are sssimed homogeneously doped , 
having acceptor concentrations and N^+. In 
addition the two regions are assumed to have different 
values of carrier lifetime and mobility. The general 
solutions for the p-base and p + region are 

Hp(s ) « A 1 cosh x + sinh x x ^ £ x £ X 

n (x) * A« cosh x + Bu sinh x X < x £ W 
P+ — * 

Application of the boundary conditions results in 
four equations and four unknowns. 

Solving for tbs unknown coefficients results in 
the following solution for the excess carriers in the 
p-base 


L -, I 

A. -2_ -J2. 1 

v D P y 



where W p + » W - X 


(l?A-i) 


As can be seen the parameter S not only describes 
the influence of tne low -high junction by the «aA*A+ 
ratio but also includes the p + region geon^try factor 
in the bracket and tte mobility and diffusion length 
associated with the two p regions. If it is assumed 
that an ohmic contact is made to the p + region, then 
the expression for S can be simplified to 





o 

p 



for a - • 


(15A-2) 


In order to determine the diode saturation current, 
the gradient of excess carriers must be determined. 
This Is obtained by differentiating equation (1UA). 
Hence 


X-X Y.v 

sinh -= — + S cosh =7— 

Mil. _ JL L p JL 

*** cosh * S sinh jj- 

V V 


(16a) 


The diffusion current is obtained by evaluating the 
gradient at x = 0; 


j - (-q) D M 

p v n dx 1 


This yields for Use p-base component of the diffusion 
current 


j ,2xi 

p V P 


S + tanh 




1 + S tanh ^ J 


Mexpev-1) 


(17A) 


or 


J pt Ip (exppy- 1 ) 


P *0 

The expression for the p-base diode saturation is 




£ + tanh 


il 


Op N.L ] V I 

A p |l + S tanh 


(ISA) 


This can be expressed in terms of -otal cell thickness 
by recognizing that W p = t-W + - (W *) where V 
and W + are the n+ diffuses region and p 4 
layer thicknesses. 


7 



Since the total diode saturation current Is com- 
posed of the n- layer and p-base contributions, an ex- 
pression for the n- layer la needed. Either of three 
expressions may be appropriate depending on the parti- 
cular model for the diffused layer. The n- layer may 
be viewed as a uniformly doped region of finite thick- 
ness W n + bounded by the space charge layer on the 
one side ar.d a metal contact of surface recombination 
velocity s on the other. For this case 





■A V 

-§-£ ♦ tanh -j- 
n n 


1 ♦ 


■A 


tanh 


(19A) 


If the n + layer Is inhomogeneous ly duped the drift 
field model should be used. Should a low high junc- 
tion model be used than an expression similar to 
equation (17A) can be derived. 


However, It will be assumed for the purpose of 
this study that the diffused region contribution Is 
always at least three orders of magnitude smaller that 
of the p-base. Simple calculations show this to be 
true assuming a homogeneous donor distribution of 
about 10 * to 10^/cc. Therefore the total diode sat- 
uration current can be approximated by (1EA) alone. 

If It is assumed that diffusion transport pro- 
cesses are dominant in a solar cell operating at open 
circuit voltage, then the V QC can be related to the 
saturation by the expression 



This equation together with (18A) is used In the study 
where = 5X10"** amp/ cm 2 and I 8C * 40 me/ cor. 


REFERENCES 

1. J. Msndelkora and J. H. Iamneck, Jr. , "Simplified 
Fabrication of Back Surface Electric Field Sili- 
con Cells and Novel Characteristics of Such Cells? 
Proceedings cf 9th Photovoltaic Specialists Con- 
ference, May 2-3, 1972. 

2. M. Wolf, "Drift Fields in Photovoltaic Solar 
Energy Converter Cells," Rroc. IEEE, vol. 51, 
p. 674, Msy 1963. 

3. J. B. Gunn, "On Carrier Accunulat ion and Proper- 
ties of Certain Semiconductor Junctions," 

J, Electron. Control , vol. 4, pp. 17-50, 1959* 

4. F. W. Dutton and H. J. Whittier, "Forward Current- 
Voltage and Switching Characteristics of p+mu 
(Epitaxial) Diodes," IEEE Trans, on Electron 
Devices, vol. ED-16, p. 458, MajT 1969^ 

3. D. P. Lieb, B. D. Jackson, and C. D. Foot, 

"Abrupt Junction Diode Theory," IRE Trans . on 
Electron Devices , vol. ED-9, P* 143, liar. T962. 

6. J. P. »!cKelvey, Solid State and Semiconductor . 

New York: Harper and Row, 1966. ’" — 

7, W. Shockley, "The Theory cf p-n Junctions in 
Semiconductors and p-n Junction Transistors," 

Bell Svst. Tech. J. , vol. 28, p. 435, 1949* 


8. D. A. Kle Inman, "Considerations on the Solar 
Cell," Bell Syat . Tech. J. , vol. 40, p. 81, 19*51. 

9. M. P. God lews ki, C. R. Baraona, and H. W. Brand - 
horst, Jr. , "Drift Field Model Applied to the 
Lithivan Containing Solar Cell," 10th Tftotovoltaic 
Specialists Conference, Nov. 13-15, 1973* 

10. Private coramnicat ion , Dr. H. W. Brandhorst , Jr., 
NASA- Lewis Research Center. 

11. F. G. Downing, J. R. Carter, Jr., and J. M. 

Denny, "The Energy Dependence of Electron Damage 
in Silicon," Proceed ings of the 4th Photovoltaic 
Specialists Conference, Vol. I, A-5-1, 19^*+* 

12. P. lies, "Effects of Processing on the Carrier 
Lifetime in Silicon Solar Cells," Contractors 
Final Report, NAS 3 -15 689, to be published. 

13. W. Zimmerman , "Measurement of Spatial Variations 
of the Carrier Lifetime in Silicon Power Devices," 
Phys. Stat. Sol, (a), vol. 12, p. 671, 1972. 

14. P. Hai-Choudhuzy, "Chemical Vapor Deposited Sili- 
con and Its Device Applications," Semiconductor 
Silicon 1973 . Edited by H. R. Huff and R. R. 
Burgess, Electrochemical Society, Inc., 1973* 

15. H. W. Brandhorst, Jr., C. R. Baraona, and C. K. 
Swartz, "Performance of Epitaxial Back Surface 
Field Cells," 10th Photovoltaic Specialists 
Conference, Nov. 13-15 » 1973* 

l£. W. van Roosbroeck, "Theory of Current -Carrier 

Transport and Photoconductivity in Semiconductors 
with Trapping." Bell Syst. Tech. J. , vol. 39. 

». 5x5-613, m/HgoT 


8 



-7G29 


EXPERIMENTAL THICKNESS-VOLTAGE PERFORMANCE OF 
BSF AND CONVENTIONAL 10 OHM-CM CELLS 



CELL THICKNESS (MICROMETERS) 
Fig. la 


EXPERIMENTAL RADIATION DAMAGE - VOLTAGE PERFORMANCE 
OF BSF AND CONVENTIONAL 10 OHM-CM CELLS 



1 MpV ELECTRON ftUENCE (e/cm 2 ) 


Fig. lb 



DEPENDENCE OF OPEN CIRCUIT VOLTAGE ON P-BASE 
DIFFUSION LENGTH FOR INFINITE BASE WIDTH MODEL 



DIFFUSION LENGTH OF p-BASE (MICROMETERS) 
Fig. 2 


DEPENDENCE OF OPEN CIRCUIT VOLTAGE ON C&L THICKNESS 
FOR RNITE WIDTH MODEL 

I$c * 40 ma/cm^, W p + = 0.5 pm, L p + = 1.0pm, p p * 10 Q-cm 
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DEPENDENCE OF OPEN CIRCUIT VOLTAGE ON ULL THICKNESS 
FOR DRIFT HELD MODEL 


INFINITE BASE 
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Fig. 4 
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EFFECT OF P + REGION DIFFUSION LENGTH ON OPEN CIRCUIT 
VOLTAGE FOR DRIFT HELD MODEL 




DEPENDENCE OF OPEN CIRCUIT VOLTAGE ON CELL THICKNESS 
FOR THE LOW-HIGH JUNCTION MODEL 
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CELL THICKNESS (MICROMETERS) 

Fig. 6 

DEPENDENCE OF OPEN CIRCUIT VOLTAGE ON NORMALIZED 
SURFACE RECOMBINATION VELOCITY S 



NG3MAI IZED SURFACE RECOMBINATION VELOCITY S 
Fig. 7 
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CALCULATED RADIATION PERFORMANCE FOR 
BSF AND CONVENTIONAL CELLS 



Fig. 8 


CROSS SECTION OF A SOLAR CELL WITH A LOW-HIGH JUNCTION 
AS PACT OF THE REAR CONTACT FOR A FiaD FREE CONDITION 


METAL CONTACT 

f777i SPACE CHARGE REGION 
n + p LOW-HIGH 



Fig. 9 
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